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INTRODUCTION 


In  recent  years,  plasticity  theory  has  been  applied  success¬ 
fully  to  a  number  of  problems  involving  soil-strength  and  soil- 
structure  interaction  (Refs.  1-5).  In  most  of  the  analyses,  a 
Mohr-Coulomb  yield  criterion  is  used  that  postulates  a  linear  re¬ 
lationship  between  the  two  principal  stresses  in  the  plane  of 
deformation.  Hill  (Ref.  6),  Sokolovskii  (Ref.  1),  and  Szymanski 
(Ref.  7)  have  set  up  the  equations  of  plastic  equilibrium  to  in¬ 
clude  general  nonlinear  yield  conditions.  Kingston  and  Spencer 
(Ref.  8),  in  applying  these  equations  to  the  problem  of  the  inden¬ 
tation  of  a  semi-infinite  mass  of  material  by  a  flat,  smooth  rigid 
punch,  have  suggested  that  the  numerical  computation  can  be  ef¬ 
fected  without  too  much  difficulty  and  with  a  relatively  small 
increase  in  computing  time. 
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SCOPE 


In  this  paper,  the  differential  equations  of  plastic  equi¬ 
librium  obtained  from  plasticity  theory  and  a  general  yield  condi¬ 
tion  are  applied  to  the  problem  of  the  interaction  between  soil 
and  a  driven  rigid  wheel.  A  computer  program,  written  in  general 
terms,  is  used  to  compute  values  of  load,  drawbar  pull,  and  torque 
for  specified  soil  parameters  and  wheel  sizes.  Results  obtained 
from  solutions,  using  a  conventional  straight  line  Mohr-Coulomb 
failure  criterion,  are  compared  to  those  in  which  a  curved  failure 
envelope  was  used  for  an  8-inch-diameter,  2-inch-wide  rigid  wheel 
being  driven  on  Jones  Beach  sand.  The  soil  parameters  for  this 
soil  were  determined  in  previous  studies  and  were  used  as  research 
tools  to  compare  theory  with  experimental  results.  This  comparison 
shows  that  the  theory  is  well  suited  for  natural  soils  having  a 
nonlinear  envelope.  Partially  saturated  soils  and  certain  types 
of  soils  with  specific  loading  histories,  e.g.,  overconsolidated 
clays,  show  failure  envelopes  that  are  significantly  curved.  In 
these  cases,  consideration  of  this  nonlinearity  is  essential  for 
the  accuracy  of  wheel-soil  interaction  prediction. 


FAILURE  CONDITIONS  IN  SOIL  BENEATH  WHEELS 


Observations  of  soil  deformations  beneath  driven  or  towed 
rigid  wheels  suggest  that  the  soil  fails  in  two  distinct  zones 
(Refs.  9,10).  The  zones  are  forward  and  to  the  rear  of  the  point 
at  which  the  maximum  normal  stress  occurs  along  the  soil-wheel 
interface.  The  reaction  of  the  soil  in  the  forward  failure  zone 
resists  the  motion  of  the  wheel,  while  the  reaction  of  the  back¬ 
ward  zone  actually  enhances  it  by  providing  a  forward  acting  com¬ 
ponent.  Under  certain  conditions,  single  zone  failure  occurs; 
however,  this  is  not  the  usual  case  and  it  will  not  be  considered 
here. 

Although  evidence  of  failure  in  the  third  dimension  exists 
(Refs.  11,12),  the  forward  and  rear  failure  zones  are  assumed  to 
form  only  in  the  radial  plane  of  the  wheel.  The  state  of  stress 
at  a  specific  point  on  the  soil-wheel  interface  is  shown  diagram- 
matically  in  Fig.  la.  The  Mohr”Coulomb  yield  criterion  was  chosen 
to  express  the  limiting  relationship  between  shear  and  normal 
stresses.  Figure  lb  shows  the  Mohr  circle  of  stress  (Circle  I) 
for  the  stress  conditions  at  the  point  on  the  interface  depicted 
in  Fig.  la.  Other  circles  (II  and  III)  can  be  constructed  for 
different  stress  states  at  other  points  along  the  interface  or  in 
the  soil  mass.  The  illustration  also  shows  how  the  major  and 
minor  principal  stresses,  their  respective  directions,  and  the 
directions  of  the  i-  and  j -characteristic  slip  lines  are  ob¬ 
tained  from  the  Mohr  diagram  for  the  stress  condition  at  the  spe¬ 
cific  point  on  the  interface.  When  these  quantities  are  trans¬ 
ferred  back  to  Fig.  la,  the  directions  of  the  principal  stresses 
relative  to  the  stresses  at  that  point  on  the  soil-wheel  interface 
are  clearly  seen. 
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Tut.  fc.Uufc.fc  envelopi  f  j  intar  or  non!  inear ;  showr  in  Fig.  lb 
.t  ct  t  ;.nuC  by  at  equation  in  which  cut  boil  strength  parameters  c 
(coLet  iouj  auc  u  (angle  of  interne  J  friction)  relate  shear  stress 
i.o  normal  fciritt  at  failure.  'J'he  nonlinear  envelope  represents 
i  he  curve  that  it  tangent  to  the  three  Mohr  stress  circles  shown. 

J  he  1 1  j  a  i.gi  i  J :  ue  usually  usee  to  approximate  the  curvilinear 
1  a  i  J.ure  envelope  <  an  introduce  errors  at  the  lower  stress  levels, 
the  fuJiure  envelope  tor  effective  stresses  is  defined  in  general 
t  .crisis  by  i  tie  equation 


in 

f  u. 
d  J 


i ' 

n 

k 


c  *  0 


(1) 


Wile  j  e 

i  -  shear  stress  at  failure 

■  j  -  effective  norma  J  stress  (as  defined  in  Ref.  13) 

n 

at  fail  ore 
i.  -  inhesion 


111,1!,  l<  -  const  ants  that  define  the  curvature  of  the 
fa i I ure  euve 1  ope . 


Km  the  lineal  case  (d  -  tit  »  I  and  k  •  cot  <k)  Eq.  (1)  re” 
dunes  to  the  more  conventional  Mohr -Coulomb  form  used  widely  in 
an  t  I  median its  (Kef.  I  i) : 


i  -  o  tan  >v  K  c  (2) 

n 

wlic re  i  -  angle  of  internal  friction  of  the  soil. 

In  the  analysis  or  the  interaction  between  a  rigid  wheel  and 
so i l ,  all  stresses  along  the  so l l -wheel  Interface  and  throughout 
the  soil  mass  are  assumed  to  act  in  planes  parallel  to  the  radial 
plane  or  the  wheel.  the  magnitude,  direction,  and  distribution  of 


stresses  beneath  a  rigid  wheel  are  functions  of  the  strength  param¬ 
eters  (c,<P,y)  of  the  supporting  soil  and  certain  properties  of 
the  wheel  [geometry,  coefficient  of  friction  (tan  6)  between 
wheel  material  and  soil,  mode  of  propulsion].  In  wheel -performance 
studies,  relationships  are  established  among  the  so-called  "per¬ 
formance  parameters"  —  load,  torque,  and  drawbar  pull  (or  drag). 
The  summation  of  the  vertical  components  of  the  normal  and  shear 
stresses  along  the  soil-wheel  interface  perimeter  and  across  the 
wheel  width  equals  the  load;  similarly,  the  summation  of  the  com¬ 
ponents  parallel  to  the  surface  is  a  measure  of  the  drawbar  pull. 
The  torque  is  determined  from  the  magnitude  and  distribution  of 
the  shear  stresses. 

The  application  of  this  method  of  wheel  performance  analysis 
is  strictly  valid  only  in  cases  where  soil  inertia  forces  can  be 
neglected  and  where  the  soil  strength  parameters  are  referred  to 
effective  stresses.  Inasmuch  as  the  kinematics  of  the  problem  are 
not  considered  in  the  analysis,  the  results  must  be  a  lower  bound 
solution.  Experiments  performed  at  Grumman  and  elsewhere  have 
shown  that  in  loose,  frictional  soils,  the  actual  failure  condi¬ 
tions  correspond  quite  closely  to  those  implied  by  the  lower  bound 
solution  (Ref.  14). 


THE  NATURE  OF  THE  FAILURE  ENVELOPE 


The  results  of  triaxial  compression  tests  on  air-dried  Jones 
Beach  sand  are  shown  in  Fig.  2.  The  origin  (0,0)  must  be  con¬ 
sidered  a  data  point  for  the  failure  envelope  of  this  material 
because  this  sand  is  not  cohesive  in  an  air-dried  condition.  Even 
without  this  constraint,  a  straight  line  cannot  be  drawn  tangent 
to  all  the  circles.  Two  limiting  straight  line  envelopes  are 
shown,  one  at  an  angle  (<p)  of  41  degrees  for  the  low  stress 
levels,  the  other  at  36  degrees  for  the  higher  stress  levels.  An 
"average"  straight  line  envelope  (not  shown)  would  tend  to  under¬ 
estimate  the  strength  of  the  material  at  low  stress  levels  and 
overestimate  it  at  high  stress  levels.  Because  it  diverges  from 
the  correct  strength  envelope  at  the  higher  stress  levels,  the 
"average"  straight  line  becomes  a  progressively  worse  approxima¬ 
tion.  A  nonlinear  envelope  corresponding  to  Eq.  (1)  is  also  shown 
in  Fig.  2.  It  is  clear  from  the  figure  that  it  provides  a  good 
approximation  to  the  Mohr  circles  at  all  stress  levels. 

The  implications  of  using  the  nonlinear  envelope  in  the  solu¬ 
tion  of  the  differential  equations  of  plasticity  and  the  associated 
stress  distributions  beneath  rigid  wheels  is  best  shown  with  ref¬ 
erence  to  Fig.  3.  In  Fig.  3a,  two  Mohr  circles,  A  and  B,  cor¬ 
responding  to  the  nonlinear  failure  envelope  are  shown  together 
with  a  cp  =  41°  envelope,  a  linear  envelope  corresponding  to 
Circle  A  for  a  cohesionless  material.  In  plasticity  theory, 
normal  stresses  are  related  to  the  intersection  of  the  failure  en¬ 
velope  with  the  a  axis.  In  the  case  of  a  cohesionless  material, 
this  intersection  coincides  with  the  origin  of  the  r  -  a  coordi¬ 
nate  system.  For  the  same  material,  the  use  of  a  nonlinear  en¬ 
velope  results  in  a  tp  intercept  on  the  a  axis  (refer  to  Fig.  3b). 
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^en  though  the  Mohr  circle  for  the  stresses  at  points  A 
and  D  in  Figs.  3b  and  3d,  respectively,  is  the  same  for  both 
the  linear  and  nonlinear  failure  envelope  (Circle  A),  the  shear 
and  normal  stresses  developed  at  the  interface,  assuming  the  same 
angle  of  interface  friction  (6),  are  different  in  the  two  cases 
(refe;r  to  points  a  and  d  in  Figs.  3a  and  3c,  respectively). 
Since  shear  stresses  at  the  interface  are  the  major  contributors 
to  drawbar  pull,  the  importance  of  using  a  nonlinear  envelope  for 
the  soil-wheel  interaction  problems  is  evident. 
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FORMULATION  FOR  A  DRIVEN  RIGID  WHEEL 


The  following  set  of  differential  equations  represents  a  theo¬ 
retically  rigorous  formulation  to  the  general  two  dimensional  prob¬ 
lem  of  determining  the  plastic  stress  states  and  slip  line  fields 
in  a  soil  mass  due  to  loading  conditions  on  specified  boundaries 
of  the  soil  mass: 


dz  =  dx  tan(0  ±  u) 

da  ±  2a  tan  9  do  -  c"0g"  "q,  js in  ( e  ±  cp)dx  +  cos(e  ±  9)dzl  =  0 


(3) 


where 

c  =  cohesion 

x,z  =  coordinates 

y  =  unit  weight  of  soil 

0  -  angle  between  x  axis  and  major  principal  stress 

9  =  angle  of  internal  friction 

i/  *  c  cot  9 

U  *  tt/4  -  <p/2 

a^  =  major  principal  stress  (see  Fig.  1) 

a^  =  minor  principal  stress  (see  Fig.  1) 

a  =  (o1  +  a3)/  2  +  ^ 

The  equations  are  derived  from  the  fundamental  equations  of  plastic 
equilibrium  and  the  Mohr-Coulomb  failure  criterion;  this  derivation 

maybe  found  in  detail  elsev*.;ere  (Refs.  15,1)-  The  coordinate  sys¬ 

tem  and  notation  for  the  specific  problem  of  a  driven  rigid  wheel 
are  shown  in  Fig.  4.  The  shaded  zones  in  the  figure  represent 
those  areas  in  which  the  shear  resistance  of  the  soil  has  been 
fully  mobilized  and  where  failure  is  occurring.  In  Eqs.  (3),  the 
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upper  sign  refers  to  the  family  of  slip  lines  corresponding  to  the 
first  characteristics  of  the  differential  equations,  the  i-lines; 
the  lower  sign  refers  to  the  second  characteristics,  the  j -lines. 
Figure  5  shows  diagrammatical ly  the  orientation  of  the  i  and  j 
characteristics  (slip  lines)  and  the  indices  used  in  computing 
values  of  x,  2,  a,  0,  and  <p  at  the  nodal  point  i,  j . 

As  no  closed  form  solution  to  these  equations  exists,  numeri¬ 
cal  methods  must  be  employed  in  which  Eqs.  (3)  are  replaced  by 
four  difference  equations.  An  iterative  scheme  is  used  to  obtain 
values  for  the  desired  variables  at  the  point  i,j  from  known 
values  either  previously  computed  at  two  adjacent  points  or  speci¬ 
fied  at  a  boundary. 
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METHOD  OF  ANALYSIS 


To  compute  the  stress  distribution  at  the  soil-wheel  interface 
and  the  slip  line  field  in  the  soil  beneath  a  driven  rigid  wheel, 
certain  soil  and  wheel  parameters  must  be  specified.  The  soil 
parameters  are:  unit  weight  y;  the  constants  n,  d,  k,  '  and  c 
for  Eq.  (1)  that  provide  the  relationship  between  shear  stress  and 
normal  stress  on  the  failure  envelope;  and  the  friction  parameter, 

5.  A  relationship  between  5  and  slip  (i)  proposed  by  Janosi 
and  Hanamoto  (Ref.  18)  is  used  in  this  study  with  the  modification 
that  a  constant  for  the  threshold  slip  (iQ)  at  which  movement 
starts,  is  included: 

-d+y/K 

tan  6  =  tan  cp  (^e  J  (4) 

where 

i  -  slip 

i  *  a  constant  to  account  for  some  threshold  perimeter 
shear  that  may  exist  before  wheel  movement  starts 

K  »  an  empirical  constant. 

The  wheel  parameters  that  must  be  specified  prior  to  the  com¬ 
putation  are:  the  wheel  radius  (rQ)  and  width  (bQ) ,  and  either 
the  entry  angle  (ag)  or  the  rear  angle  (ar)  (Fig.  4)  depending  on 
whether  the  forward  or  rear  field  is  being  computed.  For  the  analysis 
presented  here,  the  separation  angle  (am)  is  assumed  approximately 
equal  to  the  angle  6.  The  literature  contains  some  experimental  and 
theoretical  justification  for  this  assumption  (Refs.  16,17). 

Within  the  computer  program,  there  is  provision  for  definition 
of  additional  geometric  as  well  as  stress  boundary  conditions. 

For  example,  the  terrain  may  be  made  to  slope  at  an  angle  (±e) , 
and  surcharge  loadings  (w)  may  be  applied  to  the  soil  surface. 
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In  the  present  analysis,  a  level  surface  essentially  free  of  sur 
charge  is  assumed. 


The  method  of  analysis  is  based  on  a  systematic  application 
of  the  recurrence  relationships  implicit  in  Eqs.  (3)  to  the  wheel 
performance  problem  with  its  specific  boundary  condition.  In  this 
respect,  it  is  similar  to  the  method  described  in  detail  elsewhere 
for  another  problem  (Ref.  19)  with  the  added  complication  that  an 


iteration  procedure  must  be  applied  not  only  to  a 
but  to  cp. .  and  ip.,  as  well.  The  following  par, 


L.  UC  appii.CU  uut  Uili-y  LU  lj  aiiw  V  i j  , 

as  well.  The  following  paragraph  contains 


l/M  W  V>  V  T  ^  j  ^  j  u «  n  V  4.  4.  t  ilkw  4.  v/  A.  4.  V/  TV  4.11  ^  VW4I 

a  brief  description  of  the  over-all  computational  procedure. 


The  rear  field  is  computed  first  because  experience  has  shown 
that  for  given  wheel  geometry  and  soil  properties,  the  gradient  of 
the  stresses  does  not  vary  much  with  changes  in  a  .  This  is  a 
desirable  feature  in  a  trial-and-error  solution  procedure  that  re¬ 
quires  matching  the  value  of  variables  at  a  point  common  to  two 
different  computation  sequences.  The  stress  state  and  slip  line 
field  in  the  passive  zone  of  the  rear  field  are  computed  by 
Eqs.  (3)  starting  at  the  surface  for  assumed  values  of  and 

Lr  (see  Fig.  4).  In  the  radial  shear  zone,  the  same  equations 
are  used,  but  special  consideration  is  given  to  the  point  where 
the  j -lines  converge.  This  point  is  a  degenerated  slip  line, 
where  0  changes  from  the  value  at  the  passive  zone  boundary  to 
that  specified  at  the  active  zone  boundary.  The  total  change  in 
9  is  divided  by  the  number  of  slip  lines  converging  at  this  point 
to  obtain  an  equal  A 9  increment  between  two  adjacent  slip  lines. 
The  a  value  for  each  increment  is  computed  from  the  equation 
o  =  aQ  exp [2(0  -  0Q)tan  cp],  which  is  the  solution  to  Eqs.  (3)  if 
both  dx  and  dz  vanish.  The  same  equations  are  used  in  the 
active  zone  except  that  for  points  along  the  soil-wheel  interface, 
the  values  of  0  are  assigned.  Along  this  boundary,  the  values 
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of  6  and  the  wheel  curvature  must  be  considered  in  the  specifics 
tion  of  the  6  values.  The  numerical  computation  is  completed  at 
i  =  h,  j  =  k  (see  Fig.  5).  At  this  point,  the  separation  angle 
(am)  is  compared  to  6.  If  the  two  values  are  within  a  specified 
tolerance  of  each  other,  the  computation  is  continued  to  the  for¬ 
ward  zone.  Otherwise,  a  new  value  is  assumed  for  and  the 

stress  state  and  slip  line  field  computations  are  repeated  for  the 
rear  field.  This  trial-and-error  procedure  is  continued  until  the 
slip  line  field  "closes"  on  the  soil-wheel  interface  at  an  angle, 
am>  approximately  equal  to  6^  In  cases  where  closure  cannot 
be  obtained  by  varying  Lr  alone,  a  is  also  varied;  however, 
the  solution  is  not  very  sensitive  to  variations  in  a^. 

In  the  computation  of  the  forward  zone,  the  same  basic  proce¬ 
dure  is  followed:  values  of  entry  angle  (ae)  and  are  as¬ 

sumed  and  the  passive  (P) ,  radial  (R),  and  active  (A)  zone 
stress  states  and  slip  line  fields  are  computed  in  that  order  (see 
Fig.  4).  However,  an  additional  criterion  must  be  met  for  satis¬ 
factory  "closure"  of  the  forward  field.  Not  only  must  the  separa¬ 
tion  angle  am  equal  6h  k  as  obtained  from  the  rear  field  com¬ 
putation,  but  the  magnitude  and  direction  of  the  resultant  stress 
at  the  point  of  juncture  of  the  forward  and  rear  fields  must  also 
be  equal  to  that  obtained  from  the  rear  field  computations. 
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RESULTS 


All  results  relate  to  an  8-inch-diameter,  2-inch-wide,  rigid 
wheel  being  driven  over  medium  dense  (100  pcf)  ,  dry,  Jones  Beach 
sand.  The  failure  criteria  used  in  the  analyses  are  determined 
from  Fig.  2  and  are  expressed  by  Eq.  (1)  as  follows: 


JL 


1 . 1504 


JL 


1.3764 


JL 


1.1605 


-0  =  0 

(<P  =  41°) 

(5) 

-0  =  0 

(<P  *  36°) 

(6) 

-0  =  0 

(<P  -  variable) 

(7) 

where  x  and  a  are  in  units  of  pounds  per  square  foot. 

Values  of  friction  parameter  (6  *  15°,  6  =  21°,  and  6  =  25°) 
were  combined  with  the  internal  friction  angles  considered  to  yield 
a  range  of  slips  from  18  to  52  percent.  Slips  were  calculated 
from  6  and  qp  by  Eq.  (4)  with  empirically  determined  values  of 
the  constants  K  and  iQ. 

All  computations  were  performed  in  real  time  on  an  On-Line 
Systems  computer  based  on  Digital  Equipment  Corporation's  PDP-10 
processors.  After  the  computation  of  their  coordinates,  slip  line 
fields  and  stress  distributions  were  plotted  in  a  matter  of  seconds 
on  an  Adage  Inc.  visual  display  CRT  via  a  graphics  terminal  con¬ 
nected  with  the  computing  system.  This  interaction  between  the 
graphics  and  computing  systems  allowed  almost  imme.’iate  viewing  of 
the  results  to  ascertain  that  no  irregularities  existed  in  the 
solution,  and  to  acquire  a  visual  perception  of  the  interrelation¬ 
ship  among  the  many  parameters. 
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Slip  Line  Fields  and  Stress  Distributions 

Figure  6  shows  the  slip  line  fields  for  the  upper"  and  lower¬ 
boundary  straight-line  failure  envelopes  for  Jones  Beach  sand 
(cp  =  36°  and  9  -  41°)  and  for  the  nonlinear  (9  “  variable)  en¬ 
velope.  The  results  are  for  a  soil-wheel  interface  friction  of 
15  degrees.  For  comparison,  a  rear  angle  of  5  degrees  was  chosen. 
It  is  clear  from  Fig.  6  that  the  entry  angle  (a^)  in  each  case  is 
about  the  same  and  equal  to  approximately  22  degrees.  Therefore, 
the  sinkage  given  by  S  *  r  (1  -  cos  a^)  is  also  approximately 
the  same  in  each  case  (0.3  in.).  What  is  clearly  different  is  the 
geometric  extent  of  the  slip  line  field;  it  is  the  largest  for 
9  =  41  degrees  and  about  the  same  for  9  *  36  degrees  and 
9  -  variable.  This  means  that  in  the  former  case,  more  soil  mate¬ 
rial  is  affected  and  more  frictional  resistance  is  mobilized  than 
for  the  other  two  cases.  Consequently,  the  vertical  load  that  the 
wheel  can  carry  is  greater  for  9  =  41  degrees  than  for  the  other 
two  cases.  A  summary  of  the  performance  characteristics  of  the 
wheel  for  the  three  cases  is  given  in  Table  1.  It  can  be  seen 

Table  1 

SUMMARY  OF  PERFORMANCE  CHARACTERISTICS 
FOR  DRIVEN  RIGID  WHEEL  ON  JONES  BEACH  SAND 


9 

(deg) 

6 

(deg) 

a 

r 

(deg) 

Lr 

af 

(deg) 

Lf 

S 

(in.) 

W 

(lbs) 

DB 

(lbs) 

T 

(ft-lbs) 

41 

15 

5 

12 

23.1 

2 

0.32 

9.81 

0.61 

0.85 

36 

15 

5 

8.2 

22.0 

1.4 

0.29 

3.96 

0.29 

0.34 

Var. 

15 

5 

9.4 

22.0 

1.6 

0.29 

5.34 

0.41 

0.47 

from  this  table  that  although  the  extent  of  slip  line  fields  for 
cp  ■  36  degrees  and  «p  -  variable  looks  about  the  same,  the  wheel 
performance  characteristics  in  each  case  are  quite  different. 

That  the  values  of  load,  drawbar,  and  torque  are  greater  for 
CP  -  variable  than  for  cp  ■  36  degrees  is  directly  attributable 
to  the  fact  that  at  points  within  the  slip  line  field  where 
stresses  are  relatively  low,  cp  in  the  former  case  is  larger  than 
36  degrees.  Consequently,  at  those  points  more  volume  is  affected 
and  more  frictional  resistance  is  mobilized  than  would  be  the  case 
for  cp  a  36  degrees.  This  occurs  mainly  along  the  free  surface 
boundaries  and  within  the  active  zones  so  that  the  active  zones 
for  cp  -  variable  are  slightly  larger  than  those  of  cp  =  36  de¬ 
grees.  Because  the  curvature  of  the  failure  envelope  of  Jones 
Beach  sand  is  not  extreme,  this  difference  is  rather  small  and 
cannot  be  easily  noticed  in  Fig.  6.  However,  a  comparison  of  L 
values  in  Table  1  shows  that  the  effect  is  real  and  does  indeed 
exist.  From  Fig.  6,  it  is  seen  that  along  with  the  soil-wheel 
interface,  the  passive  zones  in  each  case  are  almost  identical. 

This  is  expected  since  the  stresses  within  those  zones  are  at  a 
level  where  qp  -  variable  is  almost  equal  to  36  degrees. 

Figures  7,  8,  and  9  represent  the  normal  and  shear  stress  dis¬ 
tributions  along  the  soil-wheel  interface  that  correspond  to  the 
three  slip  line  fields  shown  in  Fig.  6.  The  shape  of  the  distribu¬ 
tions  is  characteristic  of  that  determined  experimentally  for  a 
driven  wheel  by  us  and  Sela  in  Ref.  16.  The  maximum  normal  and 
shear  stresses  occur  as  specified  at  an  angle  (am)  equal  to  6  (see 
Fig.  4).  The  integration  of  the  vertical  and  horizontal  components 
of  the  normal  and  shear  stress  over  the  contact  area  yields,  re¬ 
spectively,  the  values  for  load  and  drawbar  given  in  Table  1.  The 
torque  is  determined  directly  from  the  shear  stress  distribution. 
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The  nature  of  the  lateral  distribution  has  not  yet  been  defined 
analytically  so  that  at  present  we  assume  the  normal  and  shear 
stress  distributions  to  be  uniform  over  the  width  of  the  wheel 
even  though  we  know  from  experiments  that  this  is  not  always  the 
case . 

Performance  Characteristics  As  a  Function  of  Sinkage 
and  the  Slip-Sinkage  Relationship 

The  effect  of  variable  9  on  the  performance  characteristics 
of  the  wheel  under  study  is  shown  in  Figs.  10,  11,  and  12.  In 
these  illustrations,  load,  drawbar,  and  torque  are  plotted  as  a 
function  of  sinkage  for  the  three  failure  envelopes  considered  at 
a  5  of  15  degrees.  The  illustrations  show  that  in  all  cases 
the  results  of  the  9  -  variable  analysis  lie  between  the  other 
two  and  seem  to  give  an  "average"  curve.  Closer  inspection  shows 
that  this  is  not  the  case.  At  low  sinkage  (a  situation  that  oc¬ 
curs  at  relatively  low  over-all  stress  levels),  the  9  -  variable 
curve  is  much  closer  to  the  cp  =  41  degree  curve.  This  is  ex¬ 
pected  because  at  low  stress  levels  the  cp  -  variable  strength 
envelope  approaches  the  9  ■  41  degree  Mohr-Coulomb  curve.  To 
use  an  "average"  curve  would  result  in  somewhat  lower  load  and 
drawbar,  and  higher  torque  than  predicted  by  the  9  -  variable 
analysis.  The  differences  are  small  for  the  case  in  point;  how¬ 
ever,  for  a  material  with  an  extremely  curved  failure  envelope 
the  differences  would  be  magnified,  and  use  of  either  the  9  *  41 
or  9  ■  36  degrees  failure  envelope  exclusively  would  introduce 
significant  errors  in  the  computation  of  performance  parameters. 

Figure  13  shows  the  slip-sinkage  relationship  derived  from 
Eq.  (4)  for  Jones  Beach  sand,  and  the  results  of  the  nonlinear 
failure  envelope  analysis.  From  top  to  bottom,  the  three  curves 
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shown  correspond  to  interface  friction  (6)  values  of  25°,  21°, 
and  15° ,  respectively.  Since  6  and  the  normal  and  shear 
s  tresses  acting  along  the  soil-wheel  interface  are  related  (refer 
to  Fig.  3),  it  is  clear  that  slip-sinkage  relationships  of  the 
form  shown  in  Fig.  13  are  actually  an  expression  of  soil-wheel 
interaction  and  its  effect  on  the  magnitude  and  direction  of  the 
normal  and  shear  stress  distributions. 

Relationship  of  Performance  Coefficients  to  Slip 
for  Jones  Beach  Sand 

In  the  prediction  of  tire  performance,  it  is  customary  to 
plot  experimentally  determined  performance  coefficients  versus  a 
dimensionless  number  called  the  mobility  number.  For  sand,  the 
mobility  number  is  a  function  of  the  cone  index  gradient,  vertical 
load,  and  loaded  and  unloaded  tire  geometry  (Ref.  20).  Because 
the  mobility  number  is  zero  when  there  is  no  wheel  deflection, 
this  term  cannot  be  used  as  a  predictive  parameter  for  rigid 
wheels.  Consequently,  for  this  study  on  Jones  Beach  sand,  two 
performance  coefficients,  the  pull  coefficient  (drawbar /vertical 
load)  and  the  torque  coefficient  [torque/ (vertical  load  x  wheel 
radius)]  are  plotted  in  Fig.  14  versus  another  dimensionless  param¬ 
eter,  slip.  In  the  computation  of  load,  torque,  and  drawbar,  a 
nonlinear  failure  envelope  was  used.  The  plot  differs  from  the 
conventional  pull  or  torque  versus  slip  plot  (see  e.g..  Ref.  21) 
in  that  it  contains  no  pronounced  curvature  with  increasing  slip 
that  would  correspond  to  a  maximum  pull  point  or  a  maximum  torque 
level.  However,  the  portion  of  the  plot  about  the  origin  can 
still  be  used  to  determine  the  "self-propelled"  point  (value  of 
slip  at  DB/W  -  0)  and/or  the  "towed  point"  (the  value  of  slip 
and  DB/W  at  T/WrQ  -  0)  . 
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The  results  of  two  preliminary  tests  (Test  Nos.  090101  and 
090102)  run  with  an  8-inch”diameter,  2”inch"wide  wheel  on  Jones 
Beach  sand  in  the  mobility  test  bin  at  Grumman  are  also  plotted  in 
Fig.  14.  The  points  are  seen  to  be  quite  close  to  the  theoreti- 
cally  derived  predictive  curves  for  the  same  material. 
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SUMMARY  AND  CONCLUSIONS 


A  theory  of  soil-wheel  interaction  based  on  plasticity  theory 
and  a  general  representation  of  the  Mohr  yield  criterion  has  been 
developed.  Computer  solutions  of  the  differential  equations  of 
equilibrium  are  used  to  obtain  the  slip  line  fields  and  associated 
interface  stresses  for  a  specific  wheel  geometry  and  three  differ¬ 
ent  representations  (two  linear,  one  nonlinear)  of  the  Mohr  failure 
envelope  of  the  supporting  soil  material.  Load,  drawbar  pull,  and 
torque  are  calculated  from  the  computed  interface  stresses  and  are 
compared  for  the  three  failure  envelopes.  For  each  of  these  per¬ 
formance  parameters,  the  values  for  a  given  sinkage  obtained  with 
the  nonlinear  failure  envelope  lie  between  those  obtained  with  the 
upper-  and  lower-boundary  straight-line  approximations. 

With  the  use  of  the  newly  developed  theory,  numerical  computa¬ 
tions  are  performed  that  clearly  show  the  effect  of  interface  fric¬ 
tion  on  the  relationship  between  slip  and  sinkage.  Prediction  re¬ 
lationships  for  pull  and  torque  coefficients  are  developed  as  a 
function  of  slip.  Experimental  results  are  shown  to  compare  favor¬ 
ably  with  the  analytically  derived  prediction  relationships. 

The  results  of  this  analysis  suggest  that  for  those  soils 
where  the  Mohr  failure  envelope  is  significantly  curved  (some  par¬ 
tially  saturated  soils  and  certain  types  of  overconsolidated  clays) , 
the  use  of  straight-line  approximations  or  "average"  envelopes  may 
introduce  considerable  error  in  the  computation  of  mobility  per¬ 
formance  parameters. 
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SHEAR  STRESS  (r)  psi 
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NORMAL  STRESS  (o-)  psi 


Fig.  2  Mohr  Diagram  —  Result  of  Triaxial  Tests 
on  Jones  Beach  Sand 
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A  =  ZONE  IN  WHICH  SOIL  IS 
IN  AN  ACTIVE  RANKINE 
STATE  OF  STRESS 
P  3  ZONE  IN  WHICH  SOIL  IS 
IN  A  PASSIVE  RANKINE 
STATE  OF  STRESS 
R  ■  ZONE  IN  WHICH  SOIL  IS 
IN  A  TRANSITIONAL 
STATE  OF  STRESS 
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Fig.  4  Definition  of  Problem  Geometry 
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Fig.  5  Diagram  of  Indexing  Scheme  Used  in  Numerical  Solution  of 
Equilibrium  Equations  (Arrows  along  Boundaries  Indicate 
Direction  of  Increasing  i-  and  j -Indices) 
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(a)  PHI  -36 


(b)  PHI  -41 


(c)  PHI -VAR 

Fig.  6  Variation  of  Slip  Line  Field  with  Friction 
Angle  (6  »  15°) 
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- NORMAL  STRESS 

- SHEAR  STRESS 


ANGLE  ON  SOIL-WHEEL  INTERFACE 
MEASURED  FROM  THE  VERTICAL 


Fig.  7  Interface  Stresses:  *  36* ,  6  ■  15 


ANGLE  ON  SOIL-WHEEL  INTERFACE 
MEASURED  FROM  THE  VERTICAL 


Fig.  8  Interface  Stresses:  <p  -  41°,  6 
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Fig.  9  Interface  Stresses:  cp  -  Variable,  6 


SINK  AGE  (in.) 


TORQUE  (ft  x  lbs) 


Fig.  12  Torque -Sinkage  Relationship  for  8-Inch-Diameter, 

2 -Inch-Wide  Rigid  Wheel  Driven  on  Jones  Beach  Sand 
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Fig.  14  Relationship  of  Performance  Coefficients 
to  Slip  for  Jones  Beach  Sand 


